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Abstract. We propose a model for the quantum harmonic oscillator on a discrete 
lattice which can be written in supersymmetric form, in contrast with the more 
direct discretization of the harmonic oscillator. Its ground state is easily found to be 
annihilated by the annihilation operator defined here, and its excitation spectrum is 
obtained numerically. We then define an operator whose continuum limit corresponds 
to an angular momentum, in terms of the creation-annihilation operators of our model. 
Coherent states with the correct continuum limit are also constructed. The versatility 
of the model is then used to calculate, in a simple way, the generalized position- 
dependent scattering length for a particle colliding with a single static impurity in a 
periodic potential and the exact ground state of an interacting many-body problem in 
a one-dimensional ring. 
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1. Introduction 

The quantum harmonic oscillator is a paradigmatic model with applications in all 
branches of physics, too numerous to be counted. Due to the special structure of its 
Hamiltonian, it is possible to obtain its eigenstates exactly in several different ways. 
Perhaps the most celebrated one is the algebraic solution by means of the creation- 
annihilation (ladder) operators [1], which is by far the simplest and most elegant and, 
moreover, it is the first step towards field quantization [2]. 

For more complicated problems the use of numerical methods becomes necessary. 
One of the most popular techniques is the finite-difference discretization, which is often 
employed in high-energy physics to obtain non-perturbative results [3]. However, this 
method has severe problems, since the symmetries of the original problem are usually 
lost on the lattice and can only be recovered once the continuum limit is taken which, 
in practice, is numerically impossible. As an important example, the lattice harmonic 
oscillator, though it can be written as the Mathieu equation in quasi-momentum space 
[U |5], is not factorizable and not even its ground state can be obtained in closed form. 
This fact is very disappointing, since then the supersymmetric (SUSY) structure of the 
system is not transparent - in fact, not even present - until the continuum limit is taken. 

In this article, we construct a model for the lattice harmonic oscillator which has 
a correct continuum limit. Its Hamiltonian is shape invariant [1] and, though the 
excitations cannot be accessed analytically, its ground state is exactly solvable for any 
value of the oscillator frequency and the lattice spacing. The excitations can, however, 
be obtained by solving an equation which is analogous to the Hermite equation. We 
propose then a definition of coherent states, finding that their correct continuum limit 
cannot be obtained if they are defined as eigenstates of the lattice annihilation operators, 
so their definition has to be given in terms of the displacement operator. Our model is 
completely analogous to that for a single particle in a periodic potential, and we use it to 
calculate the lowest band zero energy scattering length in a particle-impurity collision. 
We then make further use of the analogy of the model with a many-body system with 
anharmonic interactions on a finite ring which can be solved exactly for the ground 
state. 

2. Position and momentum operators on the lattice 

We define the following operators in quasi-momentum space as the momentum (p) and 
position (f ) operators. 



where d (> 0) is the lattice spacing and k G (— tt/c?, tt/c?] is the quasi-momentum. The 
operators p and x are constructed in analogy with their continuous space counterparts. 



p 



X 




(1) 



(2) 
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Note that x coincides with its continuum analog while p = hk + 0{d'^) as (i — )■ 0, and 
therefore their continuum limits are correctly described. We can write the lattice analog 
of the harmonic oscillator annihilation d and creation d^ operators as 

d ={X + iP)/V2 (3) 
d^ = {X - iP)/V2, (4) 

where the quadrature operators are defined as 

X = {mu/hf^x, (5) 
P = {mhuj)-^/^p. (6) 

However, by using the lattice operators of Eqs. ([T]) and ([2]) we see that [X, P] = i cos{kd) 
and [d, a^] = cos{kd). In other words, the canonic commutation relations are valid up 
to a factor of cos{kd). In the limit of small lattice spacing, we obtain the correct 
commutation relation of the continuous space case lim^-^-of-^, P] = i- The commutation 
relation [X, P] yields a generalized uncertainty principle (GUP) [6] of the form 

AXAP >^\{cos{kd))\, (7) 

which resembles the GUP of systems with a minimal length [7] . The GUP (j7]) does not 
imply any minimal AX (it can be zero). However, we have a maximal dispersion for 
the momentum. 



AP< AP^ax= V-^, (8) 
V mojd'^ 

which is infinite in the continuum limit, as it should. 
3. The lattice harmonic oscillator 

If we wish to construct a lattice theory for the harmonic oscillator having a similar 
structure as its continuum limit, we have to consider the operator 

N = d^d = ^{P^ + X^ + i[X,P]). (9) 

So far, the "number" operator, Eq. has exactly the same appearance as in continuous 
space. Its explicit form is given by 

The number operator written in this way looks rather unusual. If we rewrite siia^^kd) = 
(1 — cos{2kd))/2, and perform Fourier transform to direct lattice space, then we see that 
the number operator N acts as 

+ - ?(i±M+i(i^)| (11) 
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where n = x/d & Z are the lattice points. Thus, the number operator corresponds 
to a lattice with nearest-neighbor and next-nearest-neighbor hoppings with an external 
harmonic trap, plus a trivial constant. After taking the continuum limit (i — )• 0, one 
easily verifies that huN — > /2m + muj'^x^ /2 — fvjj/2. 



3.1. Ground state and lattice Hermite equation 

From now on we consider the Hamiltonian 

H = nuj{N + ]^). (12) 

If operator a annihilates a wave function in fc-space which is 27r-periodic,[8] then it is 
the ground state of H with energy Eq = huj/2. Equation d'ipo{k) = is readily solved 
and the ground state wave function has the form 

^oik) = Me-^^e^'''°<'"^\ (13) 

where ■jd = h/{mud'^). In the continuum limit, Eq. (fT3l) reduces to the well-known 
harmonic oscillator ground state, ilJoik) ~ exp(— /ifc^/mci;). We use this result to verify 
the uncertainty principle on the lattice, and find that in the ground state of the lattice 
harmonic oscillator, the uncertainty relation is also minimal. 



AXAP = - 
2 



(14) 



(V'olV'o) 

since in general AXAP > |([X,P])|/2. 

Further analogy with the harmonic oscillator in continuous space can be observed 
by solving the eigenvalue problem for the number operator Niplk) = Nil){k) with the 
ansatz il){k) = il)Q{k)(j){k). The eigenvalue problem is then transformed to the equation 

0"(fc) _ 2 A!^!lM0'(A;) + 2—Nct){k) = (15) 
moj d mtu 

that determines the unknown (f){k) for which periodic boundary conditions (PBC) 

(j){k + 2n/d) = 4>{k) are assumed. Note the analogy of Eq. fllSp with the Hermite 

equation: the naive substitution sm{kd)/d ~ k, valid for kd <^ 1, yields the well-known 

continuum limit, in which the eigenvalues N become natural numbers. 

In Fig. [T]we plot the low-energy eigenvalues of A^, Eq. (fTOj) . for a small value of the 

lattice spacing d. We see that the lowest eigenvalue is indeed zero, while the rest of the 

eigenvalues appear to be quasi-degenerate but almost linearly spaced as A'5+2 — A's = 1- 

The reason is that, in direct lattice space, the number operator includes tunneling to 

nearest and second nearest neighbors, therefore inducing the quasi-degeneracy, except 

for the ground state. The relevant eigenstates for the continuum limit are those labeled 

by even quantum numbers s and, in direct lattice space, appear to be essentially a 

superposition of the discretized Hermite functions ips/2{x = 2nd) — ips/2{x = —{2n + l)d). 
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Figure 1. Eigenvalues of N, Eq. ((TU)) . for d = — -!- 
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5.^. Coherent states 

It is natural to define now the coherent states for the lattice harmonic oscillator. First, 



we try the eigenstates of the annihilation operator, aipc 
solution of this equation is 



V'a(fc) = exp 



-i2 



h 



ak 



muo 



aipa, with a G C. The 



(16) 



with 4'o{k) being the ground state, Eq. (IT^ . If we assumed a to be any complex 
number, ipa would not fulfill the PBC Hence, if we insist on ipa to be correct, we 
have no choice but to restrict the values of a to aj = \\f^3i j G which is a 
very unsatisfactory answer since the coherent states would then be restricted to equally 
spaced real numbers. Therefore the coherent states ipa do not present a very convenient 
definition. This apparently difficult situation can be resolved in a rather elegant way, 
however, relaxing the requirement that the coherent states be eigenstates of the lattice 
annihilation operator d. To this end, we define the coherent states ^q, as solutions of 
the equation a^a(^) = a cos(/cd)\&a(/c). 



^a(^) = exp 



-i2a 



h sm.{kd) 



d 



:i7) 
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which are obviously 27r-periodic for all a G C, and have the correct continuum limit. 
We can further justify Eq. f lT7|) as a definition since even in the continuum the coherent 
states are solutions of (a + ai[X,P])'^a = 0. The only issue we cannot generahze to 
the lattice case is the usual form of the displacement operator, since on the lattice the 
Baker-Hausdorff formula is not valid due to the commutator [a, [a, a^]] 7^ 0. Therefore 
we define here the displacement (or translation) operator for lattice and continuum as 
D{a) = e~'^"^, which generates unnormalized coherent states. 

3.3. Angular momentum 

A major inconvenience of lattice discretizations, if these are introduced artificially and 
not due to a true underlying crystal structure, is the absence of continuous rotational 
symmetry. To be more concrete, we lack conservation of angular momentum or, more 
dramatically, we do not even have a definition of angular momentum on the lattice! 

We propose here a rather simple lattice analog of the angular momentum. The 
requirements this operator has to satisfy are rather relaxed: (i) it should have a correct 
continuum limit, (ii) there should be a ground state of some relevant enough Hamiltonian 
with a well defined "angular momentum" on the lattice and, (iii) the lattice angular 
momentum cannot commute with the lattice Hamiltonian. Requirement (iii) is easy to 
satisfy: take any Hamiltonian which respects the symmetries of the lattice. We discuss 
now how (i) and (ii) can be met. Consider first a two-dimensional (2D) oscillator with 
Hamiltonian = huj[Ni + iV2 + 1], with Ni = dj^aj. We define the lattice angular 

momentum in 2D, in analogy with the continuum, as 

L = X1P2 - X2P1 = i{aial - a[a2), (18) 

and we see that condition (i) is fulfilled. As promised, requirement (ii) is automatically 
satisfied by the ground state of H^'^^\ ip^^\ki,k2) = '?/'o(^i)V^o(^2), with ipoih) 
defined in Eq. ( IT3l) . and by the ground state of the 2D free particle Hamiltonian 
Hp = -2J^f^^ f,^{cos{kid) + cos{k2d)) \ki,k2){ki,k2\ , ip{ki,k2) = 6{ki)6{k2), both 
having angular momentum L = as a good quantum number. It must be noted that 
lattice angular momentum operators have been defined in the context of rotating gases 
in an optical lattice in [9], but with such definition requirement (ii) is no longer satisfied 
for the ground state of the lattice oscillator. We now show by explicit calculation that 
the angular momentum on the lattice is in general not a conserved quantity 

[L, Ni + N2] = i[ cos(/ci(i)aia| + d2d\ cos{kid) 

— did| cos(A;2(i) — cos(A;2d)a-2Q^i]; (19) 

which, as expected, is non-zero, but vanishes in the continuum limit. It must be noted 
now that the lattice angular momentum operator, Eq. f|T8|) . can be used as a definition 
not only for the model discussed here, but for any tight-binding lattice model even 
without next-nearest-neighbor hopping. 
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Figure 2. Scattering lengtlis for mV^L'^/h^ = 10~^, 1/2, 1 and 3/2 (at = from 
bottom to top). Inset: a(0) as a function of L. 



The model presented here can be applied to construct completely different systems 
and obtain some of their properties exactly. As a first application, let us consider a 
single particle moving on the real line. It is readily verified that the Hamiltonian with 
the periodic potential, 

h 



V{x) = Vosm\x/L) 



cos(x/L) 



(20) 



L V 2m 

has a ground state i^oi^) = exp[AL cos(x/L)], with = y/2mVoL/h, since the potential 
and kinetic energy operators are dual to those for the lattice Harmonic oscillator in 
quasi-momentum space. We consider a single static impurity located at xq G (— vtL, ttL] 
(this is the central site, and by translation applies to an impurity at any site), with zero 
range interaction potential Vg{x) = g6{x — xq), and we show how to get the low-energy 
scattering properties of the system in a very simple manner. First we notice that, since 
V has a purely continuous spectrum and the impurity is immobile, upon collision the 
incident waves can only acquire a phase shift. Therefore, for low-energy scattering we 
only need the periodic [ipo) and aperiodic (which we call ipj) solutions for zero energy. 
The aperiodic solution centered at the first site is given by 

^j{x) = ^jjo{x) f cia;e-'^^^°^("/^) = Vo(a;)$(a;). (21) 
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This aperiodic solution is clearly antisymmetric and it holds that = /3x + (f)p{x), 
where (f)p{x + 2ttL) = (f)p{x). Recall that without the periodic potential, this solution 
corresponds to setting 0p = 0, and the scattering length oq of a static Dirac delta 
impurity is defined |T0] by the zero-energy solution /o(a;) = 1 — |x|/ao. Clearly, /q is the 
sum of the periodic (free) solution and the aperiodic (unnormalizable) solution, with 
the inverse scattering length as a coefficient. In analogy to the free space situation, we 
define a posit ion- dependent scattering length a(xo), Xq G {—nL,nL], in terms of the 
zero-energy solution 

/(.;.„) ^^oW- l'^'W-y^'" . (22) 

a{xo) 

which, written in this way, satisfies the boundary condition 

f'i4) - f'M = 2mgf{xo)/h' (23) 
imposed by the Dirac delta, if 

»/,„<i)Vl 

(24) 



a{xo) = - 

mg 







for Xq G (— 7rL,7rL], and a{xQ + n27rL) = a{xo), with n G Z. In the simplest case of 
Xq = 0, the scattering length is shown to have the form a(0) = —h^e~'^^^/{mg). If 
xo 7^ it has to be calculated numerically. In Fig. [2]we plot the scattering length a(a;o), 
showing how strongly it depends on the position of the impurity. The scattering length 
never diverges (there is no resonance), but mga{xo)/h'^ can actually become positive 
and indeed very large with increasing Vq at Xq 7^ 0, even though the corresponding free- 
space scattering length (we assume g > 0) is negative. This means that interactions in a 
periodic potential can effectively change both quantitative and qualitatively, depending 
on where the scattering takes place. 



5. A many-body system 

As a second application, we construct a many-body Hamiltonian of interacting particles 
on a finite ring whose ground state can be obtained in closed form. We consider N 
particles on a ring of length 27rL. The position of particle i G {1, . . . ,N} is denoted 
by Xi G (— 7r/L,7rL] and its momentum by pi = —ihd/dxi, and for all the functions 
involved we use PBC We consider the following Hamiltonian 

fc2 ^ 

(25) 



2m 

1=1 

N 

dxi 



A = -^ + \L ^ sin(xij/L), (26) 

where With these definitions, the many-body Hamiltonian f l25p is 27r- 

periodic. In the limit of an infinitely long ring, L — )■ (X), we are left with N particles 
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Figure 3. Pair correlation functions p{x,x'); the clearer the color the lower its value. 
Top TV = 4 particles, bottom TV = 2 particles. Left =0.1, right XL^ = 1. Ah 
values are normalized to the peak. 



interacting via pairwise harmonic potentials. However, for any finite-size ring the 
interactions are anharmonic and contain three-body terms. 

Since the Hamiltonian H, Eq. ( I25l) . is the sum of semi-positive operators, it follows 
that H > 0. Hence, if there exists a non-singular periodic function ipQ which is 
annihilated hj all Ai, i = 1, . . . , N , then it is the ground state of H and its eigenenergy 
is zero. The set of N equations AiipQ = is easily shown to be satisfied by the wave 
function 



N 



ipoi^i, . . . ,xn) = Af Y\ ^^P 2AL^ cos(xjj/L) 



(27) 



i<j=i 
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with jV the normahzation constant. It is remarkable that for any L < oo, the ground 
state (1271) is square integrable, even if A < 0, but in taking the hmit of L — )■ oo this will 
no longer be true. 

In Fig. [3] we plot some ground state pair correlation functions, defined as 



Jn 

where = (-vr/L, Tr/L]^-^. We note that as XL'^{> 0) becomes larger, the particles 
tend to be tighter co-localized, which also happens for increasing number of particles. 

6. Conclusions 

We have constructed a lattice model of the harmonic oscillator with a correct continuum 
limit whose properties, especially in the ground state, are the perfect analogous of 
those in continuous space. We have also defined lattice coherent states and a lattice 
"angular momentum" in terms of the creation and annihilation operators of the model. 
By establishing connections with other systems, we were able to describe low-energy 
scattering in a periodic potential and an anharmonically interacting many-body system. 
These results are relevant for lattice simulations, cold collisions, many-body theory and 
quantum information [TT| [T2]. 
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